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Using a combination of first-principles theory and experiments, we provide a quantitative explana- 
tion for chemical contributions to surface-enhanced Raman spectroscopy for a well-studied organic 
molecule, benzene thiol, chemisorbed on planar Au(lll) surfaces. With density functional theory 
calculations of the static Raman tensor, we demonstrate and quantify a strong mode-dependent 
modification of benzene thiol Raman spectra by Au substrates. Raman active modes with the 
largest enhancements result from stronger contributions from Au to their electron-vibron coupling, 
as quantified through a deformation potential, a well-defined property of each vibrational mode. A 
straightforward and general analysis is introduced that allows extraction of chemical enhancement 
from experiments for specific vibrational modes; measured values are in excellent agreement with 
our calculations. 

PACS numbers: 78.30.-j, 31.15.A, 33.20. Fb, 68.43.Pq 



The ability to detect and characterize chemical species 
at the single-molecule level requires probes at the lim- 
its of present experimental resolution, and is a funda- 
mental challenge to nanoscience. Since its discovery over 
three decades ago, surface enhanced Raman spectroscopy 
(SERS) has shown significant promise for sensing individ- 
ual molecules adsorbed near metal nanostructures or sub- 
strates with nanoscale roughness [JUS]. In SERS, the con- 
version of incident light into surface plasmons near asper- 
ities on metal surfaces, combined with chemical and reso- 
nant effects, has been reported to yield Raman cross sec- 
tions increased by factors of up to 10 14 , enabling single- 
molecule detection [3HS]- While the enhancement asso- 
ciated with surface plasmons can reach 10 s [TU], remain- 
ing increases and associated changes in mode frequen- 
cies have been reported to arise from chemical adsorp- 
tion [H ITD - fll)] , as well as resonant intra- molecular and 
metal-molecule charge transfer [TTlfTi?] . Despite a wealth 
of prior theoretical studies, these chemical enhancement 
(CE) mechanisms have remained poorly understood and 
difficult to quantify [U 1121414] . While for some vibra- 
tional modes CE is estimated [H [TTJ |5D] to be about 
10-100, there is currently no clear picture for why certain 
modes are enhanced more than others. Previous theo- 
retical studies have proposed mechanisms for this mode 
dependence [Mlfl9] . but none of these models have been 
validated by experiments or more rigorous first-principles 
calculations, leaving the origin of dominant chemical con- 
tributions to SERS an open question [U IT5HT4] . 

In this Letter, we use a combination of first-principles 
calculations and experiment to demonstrate the origin of 
chemical contributions to SERS for benzene thiol (BT) 
molecules chemisorbed to Au surfaces. From our density 
functional theory (DFT) calculations of static contribu- 
tions to the Raman tensor, we elucidate the vibrational 
mode dependence [JJ IT2T - fi"4] of the chemical enhance- 



ment, explicitly relating modes with the largest CE to 
those with the greatest mode-induced shift of the molec- 
ular frontier orbital energy, as quantified through a defor- 
mation potential. Relative CE of BT vibrational modes 
measured at different probe frequencies agree quantita- 
tively with our static calculations for all binding sites 
considered. While the magnitude of CE for BT on Au is 
sensitive to binding motif, its relative value is not. 
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FIG. 1: (A) Calculated Raman spectra (from R zz compo- 
nents only) for four different binding geometries: adatom, 
hollow, adatom+H, and bridge. Inset shows our calculations 
of the orientationally-averaged gas-phase Raman spectrum of 
benzene thiol. All spectra are broadened with a Lorentzian 
of 10 cm" 1 width. (B) Binding geometries. 

DFT calculations are performed using the Vienna Ab- 
initio Simulations Package (VASP) and within a gener- 
alized gradient approximation [3TJ [55]. We model the BT 
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adsorbate-Au substrate system with an ordered mono- 
layer of one BT molecule per nm 2 bonded to a flat pe- 
riodic Au(lll) slab. Our supercell consists of 5 atomic 
layers of Au stacked along [111] with 16 atoms per layer, 
with 30 A of vacuum. The forces of the three upper layers 
and molecule are well converged to less than 1 me v/A. 
The in-plane lattice parameters are kept fixed to their 
computed Au fee bulk value of 4.17 A. A 400 eV plane- 
wave cutoff and 2x2x1 Monkhorst-Pack k-point mesh is 
used for calculations involving Au slabs. Four binding 
geometries are considered, as shown in Fig. lb: fee 
hollow (E B =0.219 eV), adatom (E B =0.446 eV), hydro- 
genated adatom (E B =0.795 eV), and bridge (E s =0.192 
eV), where Es is the calculated binding energy relative 
to a free Au surface and gas-phase C e H 5 SH in the di- 
lute limit. The gas-phase BT molecule is simulated in 
the same large supercell as the slab geometry, but using 
the r point only. 

Static Raman tensors are constructed mode-by-mode 
using a finite-differences approach, in two steps. First, 
the dynamical matrix of the system is generated by dis- 
placing each atom along each Cartesian direction by 0.03 
A. Vibrational frequencies and corresponding phonon 
eigenvectors are obtained by diagonalization of a trun- 
cated dynamical matrix treating only BT atoms, and 
Au atoms directly bonded with sulfur. (The remain- 
der of Au atoms is treated in an infinite mass approx- 
imation.) Second, we compute the static polarizability 
within a second-order finite-difference expression using a 
saw-tooth potential with a gradient of 1 mV/A, and com- 
pute its derivative as a function of the amplitude for each 
vibrational eigenmode. (See Supplementary Information 
for additional details.) Throughout the paper, all modes 
are labeled with gas-phase frequencies for simplicity. 

In Fig. la, we report Raman cross sections calculated 
from the dominant (non-evanescent) component of the 
Raman tensor, R zz , where z is the normal to the surface. 
Prominent Raman peaks in the computed spectra agree 
well with experiment [53]. Upon adsorption, BT vibra- 
tional frequencies are altered on average by 10-20 cm -1 , 
in agreement with experiments (Supplementary Informa- 
tion). From Fig. la, the presence of the Au substrate en- 
hances Raman cross sections for some modes more than 
for others. The intensities of modes with larger enhance- 
ments exhibit a stronger dependence on binding site. For 
example, while the 1100 cm -1 mode varies by more than 
two orders of magnitude across the different sites consid- 
ered, the 3151 cm -1 mode remains relatively unaltered 
and comparable to the gas-phase intensity. 

In Fig. 2a, we compare the computed charge density 
change induced by each mode, Sp(Q n ) = p{Q n ) — p(0), 
where Q n is the vibrational eigenvector. We find that 
the 1100 cm -1 mode induces significant charge redistri- 
bution within the Au surface, while the 3151 cm" 1 mode 
does not. This behavior is entirely consistent with mode- 
induced changes observed in the DFT electronic structure 
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FIG. 2: (A) Computed isosurfaces (at 1.65 x 1CT 4 e/A 3 ) of 
charge density induced by freezing-in small amplitudes (0.1 
A) of the 1100 cm and 3151 cm modes. For both modes, 
gold atoms are stationary, and charge rearrangement in the 
substrate is induced by the molecule. (B) Partial densities 
of states projected on the molecule near the Au Fermi energy 
for three different amplitudes. By inspection of the wavefunc- 
tions, the peak is the BT HOMO. 

(Fig.2b): The peak Kohn-Sham HOMO energy Ehomo 
shifts noticeably relative to the Fermi level E^ with the 
1100 cm -1 mode, whereas the 3151 cm -1 mode leaves 
the HOMO peak unchanged. Examining all modes, we 
find that modes that show larger enhancements induce a 
larger polarization response in the substrate (and larger 
shift in E f -E ff0 Mo). 

To rationalize our DFT results, we compare a model 
expression for the static polarizability and Raman ten- 
sor with our more rigorous first-principles calculations. 
We consider a single (dominant) term in an approximate 
single-particle form for the ground state electronic polar- 
izability, one that includes just a lone virtual transition 
between the HOMO, ifiHOMO , and a metallic state at 
the Au Fermi level, 4>Au- Within this two-state approx- 
imation, an adsorbate-metal interfacial contribution to 
the Raman tensor Rbt-Au can be expressed as [T5I H~4] 

Rbt A = dasT-Au w 1 dM*Mj M*M 8uj h 
0Q n UJ H dQ n u% 8Q n 

where Mi = (iphomo I Mi I 4>Au) is a dipole matrix el- 
ement in the Cartesian direction i, and u> = Epermi ~ 
Ehomo is the difference between their eigenvalues. For 
modes that are not also IR-active, terms involving deriva- 
tives of induced dipole moments can be neglected, and 

1 duo 

Eq. (1) reduces to a term proportional to ~- . 

UJ Z H dQn 

As noted previously |14j. the uijj factor in the denomina- 
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FIG. 3: (a) Partial densities of states projected on the 
molecule for different binding sites shown about the Au Fermi 



energy, (b) Difference 



R 



Adatom + H 



plotted versus 



duJH /dQ„. (c) The deviation S of these modes from linear 
trend in (b) plotted versus their contribution to the lattice 
susceptibility, which is proportional to the square of the ratio 
of the mode dynamical charge and its frequency. 



tor is consistent with a sensitivity of the overall spectral 
enhancement to the alignment of the frontier molecular 
electronic orbital with the metal Fermi energy. Indeed, 
this factor can be used to rationalize the binding site 
dependence of enhancements shown in Fig. la, and cor- 
responding partial densities of states shown on Fig. 3a. 
Calculated enhancements on the adatom site, for which 
the HOMO is close to E^, are significantly stronger than 
on other sites, where this level is broadened and further 
away from E^. However, this uijj factor is the same for 
all modes, and cannot explain the mode-dependence of 
CE computed in Fig. la. 

The strong modification of the BT Raman spectra by 
the Au substrate can be explained via computation of 
a deformation potential, Oujh / dQ n , for each mode, i.e. 
the change in molecular electronic level alignment rel- 
ative to the metallic Fermi level induced by a partic- 
ular vibration mode. In Fig. 3b, we plot the differ- 
ence R^ atom - RAdatom+H against duj H /dQ n . Because 
the Kohn-Sham HOMO level of the adatom+H geome- 
try is much further from the Fermi energy than for the 
adatom binding site (see Fig. 3a), ujh is large, and the 
adatom+H geometry has a negligible interfacial contri- 
bution. Thus, to an excellent approximation, taking the 
difference R^ atom - RAdatom+H removes intramolecu- 
lar contributions to the Raman tensor unrelated to Eq. 
(1). Indeed, in Fig. 3b, we find a remarkable correla- 
tion between the interface contribution to the deforma- 
tion potential and enhancement for most of the vibra- 
tional modes. The modes that deviate from the linear 
trend at 407 cm' 1 (Au-S stretch), 1475 cm -1 (phenyl 



ring stretch), and 3151 cm -1 (C-H stretch) have signifi- 
cant IR activity. In the static limit used here, polariza- 
tion induced by these IR active modes screens the electric 
field experienced by the molecule, leading to a reduction 
in their Raman cross sections [23] by a factor proportional 
to (Z*) /uj^, with the mode dynamical charge given by 
Z* = A/i/AQ n and where A/i is the mode-induced inter- 
facial dipole moment. In Fig. 3c, we show the deviation 
5 of the IR-active "outlier" modes from the linear trend 
observed for non-IR active modes versus (Z*) /uj 2 n . The 
correlation between 5 and (Z*) /u^ confirms that modes 
with larger contributions to the screening lead to greater 
deviations from the two-state model (Fig. 3b). We note 
however, that although important in our static calcula- 
tions, this screening effect will have impact only below 
infrared frequencies, and thus will be inconsequential for 
typical probe frequencies, where the local fields will vary 
too rapidly for the IR-active vibrations to respond. 

With this information, we can now connect the strong 
modification of Raman spectra by substrates to mode- 
specific changes of the electronic structure of the metal- 
adsorbate interface: modes with the largest interfacial 
contribution to the change in polarizability, as quanti- 
fied through a deformation potential, result in the most 
substantial chemical enhancements. For BT, these modes 
are those that break the conjugation of the HOMO (Sup- 
plementary Information). Fig. 2a shows vividly how the 
1100 cm' 1 mode breaks the resonant character of the 
carbon ring, while the 3151 cm" 1 mode leaves the ir- 
symmetry of the electrons on the phenyl ring intact. This 
suggests that for future adsorbates, the nature of modes 
with the largest CE might be intuitively rationalized a 
priori. 

To validate our theory of chemical enhancement, we 
compare with experimental Raman and SERS measure- 
ments of BT on rough Au substrates (Fig. 4a,b). 
SERS substrates, consisting of roughened SiGe surfaces 
(Fig. 4a) coated with 30 nm of Au, were incubated in 
3mM BT solution (Sigma Aldrich W361607) in methanol 
overnight, then gently rinsed with methanol and dried 
by nitrogen gas. Raman (neat solution) and SERS spec- 
tra were collected at two wavelengths, 632.8 nm and 785 
nm, using an inverted microscope set-up coupled to a 
spectrometer (Acton SpectraPro 2300i) equipped with a 
liquid-nitrogen-cooled charge-coupled device camera. To 
eliminate uncertainties associated with the number of 
molecular analytes in comparing Raman and SERS in- 
tensities, we normalize ratios of Raman and SERS spec- 
tra peak heights to the 996 cm -1 mode, which has only 
a modest enhancement from our calculations. We note 
that absolute CEs can be obtained if we normalize to 
a mode with zero deformation potential, du>H /dQ n , see 
mode 1348 cm -1 in Fig. 3b, for example. As the 1348 
cm -1 mode is not easily observed experimentally, we use 
996 cm -1 in this case. Assuming the electromagnetic 
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FIG. 4: (a) Au SERS substrate used in our measurements 
with BT. (b) Raman and SERS BT spectra, (c) Calculated 
relative enhancements versus experimental relative enhance- 
ments. Panels compare with calculations for adatom, bridge 
and hollow sites. 



enhancement is the same for all modes, this relative en- 
hancement will reflect CE. This assumption is accept- 
able for modes within a few hundred wavenumbers of the 
996 cm -1 mode, based on the width on relatively low 
Q-factor for localized plasmon resonance in Au or Ag. 

In Fig. 4c, we compare directly Raman cross sections 
from Fig. la (computed from R zz ) to the averaged so- 
lution phase data from Fig. la (inset). Binding geome- 
tries for BT on flat Au surfaces are the subject of de- 
bate in the literature; and for rough surfaces at room 
temperature, a variety of adsorption sites will be avail- 
able, and binding geometries would be subject to ther- 
mal fluctuations [9. f2"S"|. By comparing three very differ- 
ent, energy-minimized binding sites to the experiment, 
we sample different possibilities for time-averaged exper- 
imental binding morphologies. Surprisingly, all three ge- 
ometries show good correlation, with the adatom site 
showing perhaps the most linear trend. (We note that 
SERS measurements do not observe the S-H stretching 
mode at 2600 cm" 1 indicating the loss of hydrogen at 
the S-Au in the room-temperature experiments.) Impor- 
tantly, the relative theoretical enhancements, also nor- 
malized to the 996 cm -1 mode, are in excellent quanti- 
tative agreement with the present experiments, as well 
as others in literature taken from different substrates 
[2"6l |2"T] (Supplementary Information) . 

In summary, through our calculations and compari- 
son with experiment, we have demonstrated that the 
strong modification of Raman spectra by the substrate 
is a chemical effect, largely independent of laser probe 
frequency, and associated with the change in electronic 
structure of the molecule by the metal substrate. The 
mode dependence of chemical enhancement can be con- 
nected directly to interfacial contributions to the defor- 



mation potential, a well-defined intrinsic property of each 
vibrational mode and substrate. A new analysis of ex- 
perimental SERS data is introduced that allows for direct 
comparison of theoretical calculations with experimental 
data. Comparing enhancements relative to a particular 
mode, we find excellent agreement between theory and 
experiment, indicating standard DFT approaches cap- 
tures accurately and quantitatively dominant contribu- 
tions to CE, even in the static limit. Additional sup- 
port for our conclusions could be obtained by inelastic 
electron tunneling measurements, which provide direct 
access to interfacial contributions to the electron-vibron 
coupling. The quantitative connection between the de- 
formation potential and chemical enhancement provides 
new opportunities for detection and control of adsorbate- 
metal interactions through SERS. 

We thank L. Kronik, D. Prendergast, I. Tamblyn, and 
other colleagues at Molecular Foundry and UC Berkeley 
for helpful discussions. This work was supported by the 
AFOSR/DARPA Project BAA07-61 "SERS S&T Fun- 
damentals" under contract FA9550-08- 1-0257, and the 
Molecular Foundry through the Office of Science, Of- 
fice of Basic Energy Sciences, of the U.S. Department 
of Energy under Contract No. DE-AC02-05CH11231. 
Computational resources were provided by DOE (LBNL 
Lawrencium, NERSC Franklin) and DOD (HPCMP 
ARL MJM). 



[1 

[2 

[3 

K 
[5 

[6 

[7. 
[8 
[9 

[10 

[11 
[12 

[13 
[14 

[15 



Electronic address : j bneaton@lbl.gov 

M. Fleischman, P. J. Hendra, and A. McQuillan, Chem. 

Phys. Lett. 26, 163 (1974). 

D. L. Jeanmaire and R. P. V. Duyne, J. Electroanal. 
Chem. 84, 1 (1977). 

M. G. Albrecht and J. A. Creighton, J. Am. Chem. Soc. 
99, 5215 (1977). 

M. Moskovits, J. Raman Spectroscopy 36, 485 (2005). 
K. Kneipp, Y. Wang, H. Kneipp, L. T. Perelman, 
I. Itzkan, R. R. Dasari, and M. S. Feld, Phys. Rev. Lett. 
78, 1667 (1997). 

A. M. Michaels, J. Jiang, and L. Brus, J. Chem. Phys. 
B 104, 11965 (2000). 

S. Nie and S. R. Emory, Science 275, 1102 (1997). 
G. Haran, Acc. Chem. Res. 43, 1135 (2010). 

D. R. Ward, N. J. Halas, J. W. Ciszek, J. M. Tour, 
Y. Wu, P. Nordlander, and D. Natelson, Nano Lett. 8, 
919 (2008). 

K. A. Willets and R. P. V. Duyne, Annu. Rev. Phys. 
Chem. 58, 267 (2007). 

M. Moskovits, Rev. Mod. Phys. 57, 783 (1985). 

A. Campion and P. Kambhampati, Chem. Soc. Rev. 27, 

241 (1998). 

L. Jensen, C. M. Aikens, and G. C. Schatz, Chem. Soc. 
Rev. 37, 1061 (2008). 

S. M. Morton and L. Jensen, J. Am. Chem. Soc. 131, 
4090 (2009). 

E. J. Heller, R. L. Sundberg, and D. Tannor, J. Chem. 
Phys. 86, 1822 (1982). 



5 



[16] B. N. J. Persson, Chem. Phys. Lett. 82, 561 (1981). 

[17] F. J. Adrian, J. Chem. Phys. 77, 5302 (1982). 

[18] J. F. Arenas, I. L. Tocon, J. C. Otero, and J. I. Marcos, 

J. Phys. Chem. 100, 9254 (1996). 
[19] J. R. Lombardi and R. L. Birke, J. Phys. Chem. C 112, 

5605 (2008). 

[20] M. M. Maitani, D. A. A. Ohlberg, Z. Li, D. L. Allara, 

D. R. Stewart, and R. S. Williams, J. Am. Chem. Soc. 

131, 6310 (2009). 
[21] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 

Lett. 77, 3865 (1996). 
[22] G. Kresse and J. Furthmiiller, Phys. Rev. B 54, 11169 

(1995). 

[23] V. Guieu, P. Garrigue, F. Lagugne-Labarthet, L. Ser- 



vant, N. Sojic, and D. Talaga, Optics Express 17, 24030 
(2009). 

[24] M. Cardona and G. Giintherodt, "Light scattering 
in solids ii: Basic concepts and instrumentation," 
(Springer- Verlag, Berlin Heidelberg New York, 1982) 
Chap. Resonance Phenomena. 

[25] W. Zhang, B. S. Yeo, T. Schmid, and R. Zenobi, J. Phys. 
Chem. C 111, 1733 (2007). 

[26] K. B. Biggs, J. P. Camden, J. N. Anker, and R. P. V. 
Duyne, J. Chem. Phys. 113, 4581 (2009). 

[27] R. L. Aggarwal, L. W. Farrar, E. D. Diebold, and D. L. 
Polla, J. Raman Spectrosc. 40, 1331 (2009). 



